Chaotic synchronization of coupled electron-wave systems with backward waves 
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Synchronization of chaotic oscillators is one of 
the fundamental phenomena of nonlinear dynam- 
ics. In this paper we consider chaotic synchro- 
nization of two electron-wave media with inter- 
acting backward waves and cubic phase nonlin- 
earity. An efficient approach to detecting of syn- 
chronous dynamics of distributed systems is sug- 
gested based on the time scale synchronization 
concept [Chaos, 14 (3) 603 610 (2004)]. The com- 
plex space-time dynamics of the active media and 
mechanisms of the appearance of time scale syn- 
chronization are discussed. 



I. INTRODUCTION 

In last decade the investigation of chaotic synchroniza- 
tion has become a task of great interest 0, 0, 0, Hi- 
ll can be observed in a lot of chaotic oscillators of dif- 
ferent nature, including physical and biological systems 
[J IE 13 Chaotic synchronization can be used for 
secret signal transmission 0, 0, ■ It a l so seems 

very important to study the appearance of chaotic syn- 
chronization in distributed systems demonstrating space- 
time chaos Q E E E3 E EH. Four types of 
chaotic synchronization have been carried out for chaotic 
systems with lumped parameters, namely phase synchro- 
nization (PS), complete synchronization (CS), lag syn- 
chronization (LS) and generalized synchronization (GS) 
0, S 0, H3 ■ On the ranges of mentioned regimes the 
intermittence behaviour |21j takes place. To describe the 
phase synchronization the instantaneous phase 4>{t) of a 
chaotic continuous time series is usually introduced (see 
The PS means the entrainment of phases of 
chaotic signals, whereas their amplitudes remain chaotic 
and uncorrelated. Simultaneously, if Fourier spectrum of 
a signal x(t) is complex (chaotic time series are character- 



ized by Fourier spectrum without the main frequency), 
it is not always possible to describe the system dynamics 
with the help of one instantaneous phase <fr(t) (see |23ll24| 
for detail). In this case other methods can be applied to 
analyze the non-phase coherent complex systems dynam- 
ics (e.g., via curvature or via recurrence, see [25|,|2(j for 
detail) . 

In our works the new approach to the analy- 

sis of the chaotic synchronization is described, based on 
examination of continuous set of phases <p s (t) of chaotic 
time series introduced with the help of the continuous 
wavelet transform [28l I29ll30| . 

Let us consider continuous wavelet transform of time 
series x(t) 



W(s,t ) 



(1) 



where yb s ^ (t) is the wavelet-function related to the 
mother-wavelet ipo(t) as 
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In equation Q x(t) is the time series generated by a 
dynamical system. It may also contain the noise compo- 
nent. The time scale s corresponds to the width of the 
wavelet function vj sto (t) 1 and to is the shift of wavelet 
along the time axis, symbol "*" in Q denotes complex 
conjugation. It should be noted that the time scale s is 
usually used instead of the frequency / of Fourier trans- 
form and can be considered as the quantity inversed to 
it. 

The Morlet-wavelet l3ll 
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has been used as a mother-wavelet function. The choice 
of parameter value Qo = 2tt provides the relation s = 1/ f 
between the time scale s of wavelet transform and fre- 
quency / of Fourier transform. 
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The wavelet surface 

W(s,t ) = \W(s,t a )\e^ (t ^ 



(4) 



describes the system's dynamics on every time scale s at 
moment to. The value of | W(s, to)\ indicates the presence 
and intensity of the time scale s mode in the time series 
x(t) at the moment of time to. It is also possible to 
consider the quantity 



(E(s)) = / \W(s,t )\ 2 dt , 



(5) 



which is the distribution of integral energy by time scales. 

At the same time, the phase <fi s (t) — arg W(s, t) is nat- 
urally introduced for every time scale s. In other words, 
4>s(t) is a continuous function of time t and time scale s. 
It means that it is possible to describe the behavior of 
each time scale s by means of its own phase </> s (i). Let 
us consider the dynamics of two coupled oscillators. If in 
the time series x 12 (i) generated by these systems there 
is a range of time scales si < s < S2 for which the phase 
locking condition 



\4>si{t) - 4>s2(t)\ < const 



(6) 



is satisfied and the part of the wavelet spectrum energy 
in this range is not equal to zero 



«2 

E snhr = J (E{s)) ds > 0, 



(7) 



then we assert that time scale synchronization (TSS) be- 
tween oscillators takes place f24|. The condition (jfJJ) can 
be generalized to the case of m : n synchronization. In 
this case one has to use the condition 



\mcj) S ni(t) - mj) sm i{t)\ < const 



(8) 



instead of JHJl. In this case the time scales s m \ in the first 
system and s„2 in the second system would satisfy the 
relation s m 2/s n i = m/n. 

In Reference |24j the new approach was applied to the 
study of chaotic synchronization of coupled Rossler sys- 
tems and two Chua's systems. It was shown, that dif- 
ferent types of synchronous oscillations (PS, LS, ILS, CS 
and GS) may be considered as particular cases of TSS. 
It seems rather interesting to apply the concept of time 
scale synchronization to the analysis of chaotic synchro- 
nization of the distributed systems demonstrating space- 
time chaos. 

In this paper we present the results of investigation 
of chaotic synchronization in the coupled electron-wave 
systems with backward electromagnetic wave and cubic 
nonlinearity, derived with the help of the method of time 
scale synchronization. The explored system is a sim- 
ple model of microwave oscillator with backward wave 
(BWO) fin ding wi de p ractical application (see, for ex- 
ample, [32, IH, 03 HU) For similar systems the com- 
plex chaotic re gim es of space-time oscillations are typical 



The structure of the paper is as follows. Section [n] 
contains the formalism describing nonlinear interaction 
of electron beam and backward electromagnetic wave. 
Section llTTI describes the synchronization of two unidirec- 
tional coupled electron-wave systems. We demonstrate 
the efficiency of TSS method for the analysis of syn- 
chronization of spatially extended systems and discuss 
the characteristics of chaotic synchronization of such sys- 
tems. The quantitative measure of synchronization of 
spatially extended systems is described in section II VI 
In section we consider complex space-time dynamics 
of the unidirectional coupled electron- wave systems with 
the help of the concept of time-scale synchronization. 
The final conclusion is presented in section IVT1 



II. GENERAL FORMALISM 

Let us consider a simple model of two bound electron- 
wave media with backward wave and cubic phase nonlin- 
earity. 

One of the types of the active media of electron-wave 
nature is based on resonance interaction of electron wave 
and backward electromagnetic wave in the slow-wave 
structure [32, |U 03 0] . The electron wave is really 
a perturbation of density and velocity of charged parti- 
cles of electron beam propagating with the velocity vq 
along a slow- wave structure. 

In the linear approximation such system conforms to 
a model of interaction of two opposing linear waves with 
dispersion relation 



UJi = fci, 



UJ 2 



-ki 



(9) 



In this case the linear model can be formulated as a self- 
consistent system of two equations 



dF OF 



-AI, 



91 dl 

OT 0$, 



(10) 



(11) 



where an absolute instability is observed. Here F — 
\F\ exjp(jipp) and I = \I\exp(j(pj) are the normalized 
slowly variable amplitudes of interacting electromagnetic 
and electron waves, respectively, t is the normalized time, 
£ is the normalized coordinate and A is the normalized 
interaction parameter, which is proportional to the elec- 
tron beam current and the length of the active medium. 

Solution of the linear equations (|10|l and predicts 
unlimited exponential growth of amplitude of each of the 
waves. In this case it is necessary to take into account 
the nonlinear effects. One of the fundamental nonlinear 
effects is the non-isochronism of electrons-oscillators, ex- 
pressed in dependence of the frequency of the electron- 
oscillator on its energy. Non-isochronism of electrons- 
oscillators in the system of interacting waves is displayed 
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in the shift of nonlinear phase <pi of the electron wave. 
Replacing the old variables in the equations (|10|1 and Hll|) 
to new ones: r' = (t — £)/2 and £' — £ we can rewrite 
(till as 



^+3\I\^ = -A\F\e*PM<PF-<p I )), 



(12) 



where primes above new variables arc omitted. 

From the point of view of the non-isochronism effect 
the expression for dpi / dt; can be obtained from the fol- 
lowing reasons. The equation for the complete phase ipi 
of electron wave must contain the item (uj e /v e )0, defining 
the shift of the frequency of radiation due to the Dopplcr 
effect (here u e is the frequency of the electron wave and 
v e is the longitudinal velocity of its propagation). So, 
when we take into account the dependence of w e upon 
the energy of the particles W, the item proportional to 
energy W should be added at first approximation into 
the expression for a slow phase: 



fdu e \ 



(13) 



Supposing, that non-perturbed frequency uj e Q is taken 
into account in the "fast" phase, and the energy W of 
electrons (in case when their trajectories are identical) 
is defined by the energy of the wave by the expression 
W = a|/| 2 , where a is the coefficient of proportionality, 
we shall discover, that 



dpi 



da 



( duj e 



w- 



.a\I\ : 



aa 



+ V\I\ 2 , (14) 



where the derivative dip'j/da describes the phase shift of 
the wave only due to the interaction between electrons 
and electromagnetic wave F; /x = (a/v e ) (duje/dW)^ is 
the parameter characterizing the change of the electron 
velocity agreeing with the change of the electron energy 
during the interaction. We shall now define the ampli- 
tude of an electron wave as / = \I\exp(J(p'j) and deal 
with the new variable F' — F/-JJ1, I' = I / yfil, so the 
equation l|12|l will be written as follows 



-AF, 



(15) 



where the primes above new variable are omitted. 

Thus, the equation (|10fl for the backward electromag- 
netic wave will not change. 

The equations l|l(J|) and l|15fl describe the nonlinear 
oscillations in the microwave spatially extended system 
of interacting backward waves with cubic phase nonlin- 
earity. Such electron-wave system, as wtitten above, is 
a simple model of microwave oscillator with backward 
wave (BWO [H |1| and gyrotron BWO [H |H H|). 
This system is suitable for the analysis of chaotic oscil- 
lations, including chaotic synchronization, in microwave 
tubes (see, for example, BW O 1361 137L l38| and gyrotron 
backward wave generator |35l l39|). 



In an autonomous regime of BWO it is necessary to 
supplement the equations 11U|) and i|15|) with the stan- 
dard boundary conditions F(£ = l,r) = and 1(0, — 
0, r) = describing the lack of electromagnetic and elec- 
tron waves on the boundaries of the autonomous sys- 
tem. In the autonomous microwave oscillator with back- 
ward wave at7r/2 < A < 1.83 the single-frequency 
oscillations with stationary space distribution of elec- 
tromagnetic field F(a) and current 1(a) are observed. 
At A > 1.83 the multifrcquency periodic oscillations 
(reg ime of a periodic self-modulation of an output field 
|37l l38l 13^ 1 take place. At A > 2.05 the time series of 
the output field F(t,0 = 0) look like pulses with fine os- 
cillations between them, as a result of excitation of the 
complex irregular space distribution F(J[) and 1(a)- The 
reason of the latter is the fast shift of the nonlinear phase 
of the electron wave along the space coordinate of the dis- 
tributed system. At last, at A > 4.5 chaotic wide-band 
oscillations are observed. 

In present paper we consider the system of two unidi- 
rectionally coupled, spatially extended electron- wave sys- 
tems in the chaotic regime. The system is described by 
the following equations 



9Fi t 2 dF\^ 



dr 



da 



= -A X oI 



dh.2 

da 



1,2, 



-A\ ,2-Fi ,: 



(16) 



(17) 



where index "1" corresponds to the drive system and 
index "2" - to the response one. 

Unidirectional connection between microwave oscil- 
lators is brought into the system as a non-stationary 
boundary condition for the slowly varying amplitude of 
the electromagnetic field F2 of the response system thus 
the boundary condition for the first conducting fissile 
medium remains constant 



h(a = 0,r) = 0, F 1 (a = l,r)=0, 



(18) 



7 2 (£ = 0,t) = 0, F 2 (^ = 1,^=^(^ = 0^). (19) 

Hence p — Rexp[j9] is the complex coupling coefficient 
(R is the coupling strength and 9 is the phase of coupling 
coefficient). 

For drive and response systems the values of control 
parameters are taken as A\ — 4.2 and A2 = 4.9, that 
corresponds to the regimes of chaotic oscillations in the 
active electron- wave systems with backward waves. The 
highest Lyapunov exponents are A^i = 0.224 ± 0.010 
and Xa2 = 0.617 ± 0.006 for the first and second control 
parameter values, respectively. The value of the phase of 
coupling coefficient 9 does not influence the processes in 
the coupled systems and is fixed 9 — w. In present paper 
the establishment of the chaotic synchronization regime 
in two spatially extended electron-wave systems (BWO) 
is analyzed at the variation of coupling strength R. 
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FIG. 1: Power spectra and phase pictures of output electro- 
magnetic field oscillations, (a) Drive oscillator with the value 
of dimensionless parameter A = 4.2. (b) Response oscillator 
without external signal (A — 4.9 and R — 0). (c,d) Response 
oscillator with the external signal. The value of coupling pa- 
rameter R is (c) R = 0.2 and (d) R = 0.5 



III. CHAOTIC SYNCHRONIZATION OF 
SPATIALLY EXTENDED ELECTRON-WAVE 
SYSTEMS 

Let us consider the behavior of the coupled systems 
with the increase of amplitude of coupling coefficient R 
and the fixed values of control parameters (the values 
mentioned above are taken). We shall investigate the 
oscillations F r i :2 (T) — Re [-Fi,2(£ = 0,r)] at outputs £ = 
of both coupled systems. 

Figures ^ andQJ) show the power spectra and phase 
pictures (reconstructed by delay-coordinate embedding 
method ^[) of oscillations at the outputs of drive 
and response systems in autonomous regimes (i.e. at 
R = 0). 

The output field dynamics F r \ of the drive oscillator 
(which acts upon the input of the interaction space £ = 1 
of response oscillator) is chaotic (see Fig.^). In the 
power spectrum the high noise background at the level 
—30 dB is observed. There are two main spectral com- 




0.96 



0.3 0.4 0.5 " 

FIG. 2: The dependence of the synchronized time scale range 
[si;«2] on the coupling strength R for base time scales Sb 
(bottom) and s m (top) 



ponents with frequencies fb and f m in this spectrum. 
High frequency f b is close to frequency i1>/2tt of synchro- 
nism of the linear noninteracting electron and backward 
electromagnetic waves. At the decrease of dimensionless 
parameter of the electron beam A\ the regime of periodic 
oscillations with the frequency close to fb is established. 
The second frequency which is marked as f m and is set- 
tling down in a low-frequency part of a power spectrum, 
determines the low-frequency modulation of complex am- 
plitude |-F\(t)| of the electromagnetic backward wave. 

With the growth of control parameter A the oscilla- 
tions in the electron- wave system complicate, that is ob- 
viously seen from the analysis of the power spectrum and 
the chaotic attractor in Fig.^b), of the second system 
at A — 4.9 and R = (autonomous oscillations). It 
is shown in Fig.^b), that the noise background of the 
power spectrum rises up to a level —(5 -j- 10) dB. With 
such noise background it is already impossible to select 
the main spectral components. 

Let us consider the behavior of the response system 
under the influence of the drive system signal, at the 
increase of coupling coefficient R. The growth of the 
coupling strength parameter R means that the power of 
external signal acting on the response system, increases 
as R 2 . 

With the growth of a coupling parameter R the power 
spectrum of generation of the response system changes 
that is illustrated by Fig.Q]c,d). At the small R < 0.3 
[Fig.QJc)] upon the noise base of the response system 
spectrum there is a weakly delineated frequency (marked 
by an arrow), corresponding to the base frequency fb of 
the driving system spectrum [see Fig.QJa)]. At the large 
values R > 0.3 the energy of the spectral component 
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corresponding to the base frequency (marked by arrow 
in Fig.^d); R — 0.5) increases. 

It is necessary to mark, that both systems have attrac- 
tors with ill-defined phase. For these chaotic signals it is 
impossible to introduce the instantaneous phase <f)(t) of 
chaotic signal correctly. It is clear that in such case the 
traditional methods of the phase synchronization detect- 
ing (see Ref. [J, y, 112, l23|| ) fail. On the contrary, the 
approach based on introduction of instantaneous phase 
4>s if) by the wavelet transform p^L |27j can easily help 
to detect the time scale synchronization between chaotic 
oscillators. 

The dependence of synchronized time scale range 
[s\; S2\ on coupling parameter R is shown in Fig. EI The 
appearance of the synchronized time scale range corre- 
sponds to the TSS regime. 

The peculiarity of dynamics of the system of two unidi- 
rectional coupled spatially extended systems is that the 
TSS of chaotic oscillations of each of subsystems at a 
large value of coupling parameter R is observed for two 
ranges near to base time scales s& — 1/ fb and s m = l/s m , 
respectively. In Fig. [21 the ranges of phase locked scales 
are shown for the areas near to the base time scale st> 
(lower figure), and the time scale s m (upper figure). 

From Fig. Hit can be seen, that for coupling coefficient 
values R < 0.29 the establishment of a chaotic synchro- 
nization regime is not observed. The latter means, that 
there are no time scales s which would satisfy the condi- 
tions of TSS © and 0. 

Figure [3] illustrates the dynamics of phase difference 
4>si(t) — 4>s2(t) for different time scales for R = 0.2. 
It is clear, that there is no time scale range [si;s2] for 
which the phase difference 4> s i(t) — 4>s2(t) (s S [s±; S2]) is 
bounded. 

For the coupling parameter R exceeding the critical 
value R c ~ 0.29 we observe the time scale range As — 
S2 — si for which the requirement of capture of phases 
© is satisfied and the energy relating on this range 
is nonzero. The latter means, that in the coupled active 
media the TSS regime is established |2ll27l|. 

The behavior of the phase difference 4> s i(t) — 4> S 2(t) 
for this case is presented in Fig.^b). So, we can say, 
that the time scales s = 1.055 (and close to them) of 
chaotic oscillations of electromagnetic field amplitudes 
on the output of coupled systems are synchronized with 
each other. At the same time the other time scales (e.g., 
,s = 1.5) remain uncorrelated. For such time scales the 
phase locking is not observed [see Fig.^b)]. 

Let us note, that with the growth of coupling param- 
eter the TSS appears, at first, near the base scale sj, at 
R = Rb (it is marked in Fig.[5J). With the further increase 
of coupling coefficient, at R = R m > Rb the phase lock- 
ing of time scales occurs in the area of a time scale s m , 
corresponding to modulating frequency f m of amplitude 
of the output field of the response BWO. 

Let us consider at first, why the range As of synchro- 
nized time scales (i.e. time scales, for which the con- 
ditions of synchronization © and Q are fulfilled) ap- 



150 -, 



100 - 




FIG. 3: The phase difference </!>si(£) — 4>s2(i) for two unidi- 
rectionally coupled spatially extended electron-wave systems. 

(a) R — 0.2. There is no synchronization between systems. 

(b) R = 0.3. The time scales s — 1.055 (red line) of both 
systems are correlated with each other and the time scale 
synchronization regime can be observed 



pears near the time scale Sb corresponding to the base 
frequency fb = 1/sb- For this purpose we shall study 
wavelet power spectra (£^2(5)), defined from the equa- 
tion (J5J, for drive and response systems. Wavelet spectra 
of the output field oscillations of drive (a stroke line 1) 
and response (a solid line 2) BWO are presented in Fig.0] 
for different values of coupling parameter R. Let us mark, 
that the wavelet power spectrum (Ei(s)) of the drive 
system does not vary (since the value of controlling pa- 
rameter Ai is fixed and the system is autonomous) and 
is presented in all the figures, for different values of R, 
as a matter of convenience comparison with the wavelet 
power spectra (E2) of the response system. 

In Fig.QJa) wavelet power spectra of autonomous os- 
cillations are shown (R = 0). One can see that the os- 
cillations in each of the spatially extended systems have 
essentially different wavelet power spectra. So, in wavelet 
spectrum (Ei(s)) of the drive system the base time scales 
Sb (corresponding to the high frequency /),) and s m 
(corresponding to modulation oscillations with the fre- 
quency /,„ = l/s m ) stand out against the background. 
In wavelet spectrum of an autonomous response system 
the time scales corresponding to high-frequency spectral 
component (they are marked by an arrow in Fig.^Ja)) 
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FIG. 4: The energy distribution in wavelet spectra (E(s)) 
for output field oscillations of drive (stroke line) and response 
(solid line) spatially extended systems, (a) R = 0. There are 
autonomous oscillations in the systems, (b) R = 0.2. There is 
no time scale synchronization between systems, (c) R — 0.3. 
There is time scale synchronization near the base time scale 
Sb (Synchronized scale region is marked by gray box S). (d) 
R — 0.5. Time scale synchronization near both base time 
scale Sb and modulation time scale s m is observed. 



also precisely stand out against background. However, 
the arrangement of a peak of Sb in wavelet power spec- 
trum (E 2 (s)) does not correspond to a base time scale Sb 
of the drive system. As to large time scales (correspond- 
ing to low modulation frequencies), their wavelet power 
spectra are essentially different. The wavelet spectrum 
of the response system in the area of the time scale s m is 
continuous and homogeneous, i.e. the modulation oscilla- 
tions are characterized by continuous noise-like spectrum. 
For small values of coupling coefficient R, the shape of 
wavelet power spectrum changes [see Fig.0fb)). How- 
ever, time scale synchronization is not observed, which is 
apparent from Fig.^a). 

With the increase of R the conditions of TSS JBJ and 
@ are satisfied only by the time scales s whose energy in 
power spectra of the drive and response systems is rather 
large. Really, at the values of coupling coefficient R » 
0.3 as it was already considered above, the conditions of 
TSS are fulfilled only for time scales near the basic scale 
Sb = 1/fb ~ 1-0 (see also Fig.[2J|. As it follows from 
Fig-EIc), in the range of time scales marked in Fig.^Jc) 
by grey color and symbol S, the energy of wavelet power 
spectrum is maximum. 

From Fig.^Jc) one can also see that the energy of 
wavelet spectrum, being fallen on the time scale s m of 
the response system simultaneously increases. However, 
at R w 0.3 its energy is not large enough, and TSS of 
field modulation oscillations is not observed. 

At large coupling coefficients R > 0.35 the shape of 
wavelet power spectra of output field oscillations of re- 
sponse and drive systems become close to each other (see 
FigJUd)). Base time scales s m and Sb are presented in 
the power spectra with approximately equal intensity. 
The analysis of time scales dynamics shows, that the con- 
ditions of synchronization @ and J7J are fulfilled for both 
most intensive time scales in wavelet power spectra and 
two ranges of time scales, in which phase locking occurs 
[see Fig.EJfor R > 0.35], are accordingly selected. 

Thus, the TSS of two coupled chaotic spatially ex- 
tended systems is observed first of all for those time scales 
s whose energy in the chaotic power spectrum is great. 

With the increase of coupling parameter the range of 
synchronized time scales [s± , S2] , for which the TSS condi- 
tions 10 and (J7J are satisfied, is expanded. Time scales 
which are close to the most intensive scales in wavelet 
power spectrum and whose energy is large enough, start 
to be involved in synchronous dynamics. 

The latter is well illustrated by Fig. [21 from which it 
can be seen, that for both base time scales s& and s m 
extension of the range of phase locked scales As with the 
increase of R is observed. At large values of coupling coef- 
ficient R the range of synchronous scales As& = S26 — Si& 
in the area of the base scale s& ceases to be expanded 
and, moreover, is narrowed down at R > 0.5. At the 
same time in the area of scales s m , describing modula- 
tion of amplitude of the output field of BWO, the range 
of synchronous scales As m = s 2 m — Sim increases lin- 
early with the growth of i?, at the expense of growth of 
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top bound of area of synchronized scales S2m- At last 
we must note that with growth of R, the energy of the 
wavelet power spectrum of a response system is redis- 
tributed in such a manner that the increase of energy of 
larger scales s > s rn is observed. Then the most inten- 
sive time scales in wavelet spectrum appear locally, near 
the top bound si m of synchronized time scales, therefore 
these time scales are first to be synchronized with the 
increase of coupling coefficient R. 



IV. MEASURE OF SYNCHRONIZATION OF 
SPATIALLY EXTENDED ELECTRON-WAVE 
SYSTEM 

Introducing continuous set of time scales s and the 
instantaneous phases associated with them, as well as 
separation of the range of synchronous time scales As = 
s 2 — S\ allows us to inject the quantitative performance of 
a measure of chaotic synchronization of coupled systems 
pil l^ . This measure 7 can be defined as a part of 
wavelet spectrum energy being fallen on the synchronized 
time scales 

J=^J(E 2 (s))ds, (20) 

Sin 

where [s m ;s{,] is the range of time scales for which the 
condition JBJ is fulfilled and E 2 is the full energy of the 
wavelet spectrum 

+00 

E 2 = f(E 2 {s))ds (21) 


of the response system. 

The energy distribution over different time scales is a 
very useful characteristics for complex system behavior 
description. It is widely used in different applications 
(see e.g., [HIHI). 

This measure 7 is equal to zero for the nonsynchronized 
oscillations. 7^0 means that in the coupled systems 
the conditions of TSS © and J7J) are implemented. The 
value 7 = 1 shows, that oscillations in each of subsystems 
are identical or shifted upon some lag of time T T . This 
is a complete or lag chaotic synchronization £| . Increase 
of 7 from up to 1 testifies to the increase of a part of 
the wavelet spectrum energy, being fallen on synchronous 
time scales s. Actually, the measure 7 characterizes how 
close to each other the chaotic oscillations are in both 
spatially extended systems. 

In our case when the synchronous behaviour is ob- 
served for two base scales s m and Sb, it is possible to in- 
troduce measures of synchronization 7^ and 7 m according 
to each of time scales Sb and s m , and also the integral per- 
formance of the TSS measure of the systems: 7 = jt+Jm- 
Dependence of the measure 7 on the coupling coefficient 




0.25 0.35 0.45 0.55 R 



FIG. 5: The dependence of the synchronization measures 7f,, 
7 m and 7 of the response system on the coupling strength R 

R characterizes degree of chaotic synchronization of os- 
cillations in the electron-wave systems with a backward 
wave. 

In Fig.Elthe dependence of the synchronization mea- 
sures 7b, 7 m and 7 for the response systems upon the 
coupling parameter R is shown. It is obvious, that at 
small R < 0.4-^0.5 the main part of energy of the wavelet 
power spectrum, being fallen at synchronized time scales, 
corresponds to the base scale Sb (jb > 7m)- 

Thus, in the coupled systems there is the capture of 
phases of time scales s ~ Sb close to the time scale cor- 
responding to the frequency of synchronism of electronic 
and electromagnetic waves. With further growth of R 
the capture of frequencies close to the modulating fre- 
quency f m of the output signal of the driving BWO takes 
place. With the increase of R the range of the synchro- 
nized scales describing low-frequency chaotic modulation 
of the output signal also increases. Hence, j m (the nor- 
malized part of energy of synchronized scales correspond- 
ing to low-frequency modulation oscillations) grows with 
the increase of R. At R > 0.5 part of energy of syn- 
chronous low-frequency modulation time scales "f m be- 
comes more than the value 7^, describing part of energy 
of synchronous high-frequency oscillations. 

The integral part of energy 7, being fallen on syn- 
chronous time scales, monotonically increases with the 
growth of value of coupling coefficient R [see Fig.|3]. 
However, even at large values of R, when the response 
active medium is hardly influenced by the signal of the 
driving system (R 2 ~ 0.3 -r 0.4), the part of the energy, 
being fallen upon synchronous time scales, does not ex- 
ceed value 7 ~ 0.3. 

Thus, at chaotic synchronization of two unidirectional 
coupled BWO at the increase of coupling coefficient, we 
can observe TSS near the most intensive time scale Sb cor- 
responding to frequency fb of synchronism of electronic 
and electromagnetic waves. With the growth of R there is 
an increase of energy of synchronized scales in the area of 
the base scale Sb- At the same time, a range of synchro- 
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nized time scales As m , corresponding to low- frequency 
modulation of oscillations of field F, appears. 

With the further increase of R the growth of wavelet 
power spectrum energy, being fallen on synchronous time 
scales, is due only to the increase of the range of synchro- 
nized time scales As m . The given behaviour of two unidi- 
rectional coupled systems is quantitatively described by 
the measure of chaotic synchronization 7 = "f{R). 



V. SPACE-TIME DYNAMICS OF RESPONSE 
CHAOTIC ELECTRON-WAVE SYSTEM 

Let us consider the physical processes in the non- 
autonomous (response) active medium when the chaotic 
synchronization appears. First of all, we are interested 
in space dynamics of response active medium. 

The external chaotic signal generated by the drive dis- 
tributed system influences on the point with coordinate 
£ = L of the response active medium and then spreads 
towards electron wave to the output £ = border. The 
dimensionless length of the active media denoted as L 
is equal to 1. The synchronization of the whole active 
electron-wave system discussed above means the time 
scale synchronization of external signal F ex ti(r) and out- 
put signal F out 2{T) — F^{t, £ = 0) of response media 
takes place. To consider the space-time dynamics of re- 
sponse system we have analyzed the presence of time 
scale synchronization between the external signal and 
the time series F2 (r, £) obtained from different points 
< £ < L of interaction space. 

Let us consider the synchronization on time scales 
which are close to the most intensive scale s;, of exter- 
nal chaotic signal wavelet spectrum. The synchronization 
of these time scales appears when the value of coupling 
strength R is small enough and when the energy of exter- 
nal signal P ext — R 2 \F ext i | 2 is also small (see Section ITTll 
for detail). The range As;, of synchronized time scales 
may be considered as a function of interaction space coor- 
dinate of active medium. This function Asb — As{,(£) as 
well as the synchronization measure 7 = 7(£) (|20[1 are the 
characteristics describing the space dynamics of chaotic 
synchronization appearance when the external chaotic 
signal acts on one of the border of the distributed system. 

Fig-EIa) shows areas of synchronized time scales on the 
plane "coordinate (L — £) — time scales s" for three dif- 
ferent values of coupling strength R. In Fig. Eta) one can 
see that the width of the synchronized time scales area 
is maximum in the points of interaction space which are 
close to the border of system on which the chaotic signal 
of the drive active medium acts (see the border condi- 
tions (|19|l for the field F2 of the response system) . When 
the coordinate £ decreases, the area of synchronized time 
scales As;,(£) converges gradually. The measure of syn- 
chronization 7i)(£) behaves along the interaction space in 
the similar way [see Fig. HJb)] . With the decreasing of £- 
coordinate the part of wavelet spectra energy being fallen 
on the synchronized time scales Asb(£) also decreases. 



s 't> ~L 




FIG. 6: (a) The areas of synchronized time scales on the plane 
"coordinate (L — £) — time scales s" . The coordinate (L — f) 
is counted from the end of the response system, (b) The 
dependence of synchronization measure 75 on the coordinate 
of interaction space (L — £). The black squares (■) correspond 
to the value of coupling parameter R — 0.4, circles (o) - 
R = 0.3 and black triangles (A) — R = 0.2 

When the coupling strength R is rather small (R < 
0.29) and the synchronization of two coupled active me- 
dia is not observed [see Sections IIIII and llV) . the width 
of the synchronized time scales area becomes equal to 
zero at the point of interaction space £ s < L. It means 
that the intricate space dynamics of phase locking takes 
place in the asynchronous regime. The dependence of the 
width of the synchronized time scales area on the coor- 
dinate £ (see also Fig. Eta), case R = 0.2) is described 
by 

As ± i > £ s , , , 

As = £ < [ZZ) 

So, in this case it is possible to say that the interaction 
space of the response active medium is separated to two 
regions. 
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First of them is characterized by length L s = L — £ s 
and abuts on the border £ = L of the response active 
system. The chaotic oscillations -FifV,^ (£ € (£ S ,.L)) in 
this region are synchronous (in terms of time scale syn- 
chronization) with the external chaotic signal F ex ti(r) 
generated by the drive system. Therefore, the synchro- 
nization measure is not equal to zero: 76|c>£ 7^ 0. This 
region of interaction space has naturally been called as 
"the region of synchronous oscillation" . Its length L s has 
been called "length of synchronization" . 

In the second region £ € (0, £ s ) the chaotic synchro- 
nization is destroyed and the range of the synchronized 
time scales <|22[) as well as the synchronization measure 
is equal to zero. 

With the increasing of R the synchronization length 
also increases. For some value R = Rb the length of 
synchronization becomes equal to the length of the in- 
teraction space L s = L (see Section IIIII and Fig. [2). 
This means that in all space of response of electron- 
wave medium the chaotic synchronization takes place 
and, therefore, one can say about synchronization of con- 
sidered systems as whole. This regime of chaotic synchro- 
nization has been described above in sections IlIII and llVl 
when the time scale dynamics of chaotic signal get from 
output £ = of response system has been considered. 

Fig. HJa) corresponding to the case of synchroniza- 
tion of the whole response active medium shows the bor- 
ders of synchronized time scales s along the interaction 
space. The value of coupling strength has been selected 
as R = 0.3 w Rb- There is the range As of synchronized 
time scales for every point < £ < L. With coupling 
strength increasing the range of synchronized time scales 
also increases [see Fig.^a), R = 0.4]. This result agrees 
well with Sec. IIIII The value of wavelet power spectrum 
energy being fallen on synchronized time scales is equal 
approximately to 30 4- 40% in the point with coordinate 
£ = 0.99L and decreases to 1 5% at the output border 
£ = of the response active media. 

In Fig. [7| the dependence of the length of synchroniza- 
tion L s on coupling parameter R for the base time scale 
Sb is shown (solid line). It is clear, that with the in- 
creasing of coupling strength the length of synchroniza- 
tion grows monotonically according to the linear law: 
L S (R) ~ R. It follows, that with the growth of cou- 
pling coefficient R we have the linear increase of the re- 
gion where time scale synchronization near base scale Sb 
is observed. Such behaviour of the response system is 
determined by the features of interaction of electronic 
and electromagnetic waves. Destruction of time scale 
synchronization is defined by violation of phase ratio be- 
tween backward waves of electron current I(t, £) and elec- 
tromagnetic field F(t, £). 

So, at the fixed length of the response system, the 
phase shift Aip between field and current waves, origi- 
nating due to the influence of external signal upon the 
response system, reduces on the length of synchroniza- 
tion L s to destruction of time scale synchronization on 
some time scale s. The phase shift is defined as Atp = 




FIG. 7: Length of synchronization of response spatially ex- 
tended system as function of coupling strength R. The de- 
pendence for base time scale Sb is pointed by circles and ap- 
proximation L s ~ R is plotted by solid line. The numerical 
dependence and approximation for modulating time scale s m 
is pointed by squares and dashed line 

2tt|/& - fs\£,s/v g , where f b = l/s b , f s = 1/s and v g is 
the the group velocity of the backward electromagnetic 
wave. 

On the other hand, dependence of the phase shift Aip 
upon the amplitude of the external field R\F2\ can be 
presented as Aip = KRF2, where k is the coefficient of 
proportionality depending on length of the response spa- 
tially extended system A2 . Comparing the estimated ex- 
pressions for the phase shift Atp it is found out that the 
length of synchronization L s is proportional to R, as the 
numerical simulation shows (see Fig. 01 . 

It is important to note that curves shown in Fig. 
correspond to the base time scale s& and time scales 
close to it. The similar results have been obtained for 
the space-time dynamics of scales which are close to the 
modulation time scale s m [see Fig. [7] dashed line]. The 
difference consists in that in this case the length of syn- 
chronized time scales area (length of synchronization) be- 
comes equal to the system length for the bigger values of 
coupling strength than for the base time scale Sb- This 
also agrees well with the results discussed above. 



VI. CONCLUSION 

In conclusion, in our article we have considered 
the chaotic oscillations in two unidirectionally coupled 
electron-wave active media. The approach of chaotic 
synchronization analysis proposed in our work p4l |2^ | 
has been used. This approach is based on the consider- 
ation of continuous set of time scales and chaotic signal 
phases associated with them. The application of this 
method allows us to investigate the chaotic synchroniza- 
tion of distributed active media demonstrating the intri- 
cate dynamics whereas the traditional methods of chaotic 
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signal phase introduction fail due to the complex topol- 
ogy of the attractor. At the same time we succeed in di- 
viding the chaotic synchronization into synchronization 
of low-frequency modulation oscillations of electromag- 
netic field and high-frequency generation one. One can 
see that it was impossible to do it using traditional ap- 
proaches of chaotic synchronization analysis. So, the re- 
sults described in our work allow to understand more 
clearly the common mechanisms of chaotic synchroniza- 
tion and justify the efficiency of new TSS-approach to 
chaotic synchronization analysis. 

The main results are the following. With the increase 
of coupling strength the chaotic synchronization regime 
appears and is characterized by phase relationship be- 
tween time scales of chaotic oscillations of drive and re- 
sponse active media. First of all, the synchronization 
of time scales corresponding to high-frequency spectral 
component appears. The synchronization of time scales 
corresponding to low-frequency oscillations takes place 
with the larger values of coupling strength R. When the 
chaotic synchronization of the whole active medium is not 
observed the interaction space may be divided into two 
regions corresponding to synchronous and asynchronous 
oscillations, respectively. The length of synchronous os- 
cillation area L s grows when the coupling strength R 



increases. The appearance of chaotic synchronization of 
two unidirectionally coupled active media is determined 
by condition of equality of synchronization length and 
length of the system. 

It should also be noted that the method of investiga- 
tion of chaotic synchronization of two coupled distributed 
systems discussed above can be easily applied to the anal- 
ysis of higher spatial dimensions systems, such as 2D and 
3D ones. In this case one have to analyze time series 
x(t, r) obtained from different points r of space in the 
manner discussed above. 



VII. ACKNOWLEDGMENTS 

We express our appreciation to Professor Dmir- 
tiy I. Trubetskov for valuable discussions. We also thank 
Dr. S.V. Ercmina for language support. This work has 
been supported by U.S. Civilian Research & Develop- 
ment Foundation for the Independent States of the For- 
mer Soviet Union (CRDF, grant REC-006) and Russian 
Foundation for Basic Research (grant 05-02-16273). We 
also thank "Dynasty" Foundation. 



[1] Rosenblum M.G., Pikovsky A.S., Kurths J., Phys. Rev. 

Lett. 76, 1804 (1996). 
[2] Pikovsky A., Rosenblum M., Kurths J., Int. J. Bifurca- 
tion and Chaos 10, 2291 (2000). 
[3] Boccaletti S., Kurths J., Osipov G., Valladares D.L., 

Zhou C, Physics Reports 366, 1 (2002). 
[4] Pikovsky A., Rosenblum M., Kurths J., Synhronization: 

a universal concept in nonlinear sciences (Cambridge 

University Press, 2001). 
[5] Ticos CM., Rosa E., Pardo W.B., Walkenstein J.A., 

Monti M., Phys. Rev. Lett. 85, 2929 (2000). 
[6] Anishchenko V.S., Balanov A.G., Janson N.B., Igo- 

sheva N.B., Bordyugov G.V., Int. J. Bifurcation and 

Chaos 10, 2339 (2000). 
[7] Tass PA. et al., Phys. Rev. Lett. 90, 088101 (2003). 
[8] Prokhorov M. at al., Phys. Rev. E 68, 041913 (2003). 
[9] Murali K., Lakshmanan M., Phys. Rev. E 48, R1624 

(1994) . 

[10] Yang T., Wu C.W. and Chua L.O., IEEE Trans. Circuits 

and Syst. 44, 469 (1997). 
[11] Liu Yaowen, Ge Guangming, Zhao Hong, Wang Yinghai, 

Phys. Rev. E 62, 7898 (2000). 
[12] Garcia-Ojalvo J., Roy R., Phys.Rev.Lett. 86, 5204 

(2001). 

[13] Kurths J., Pikovsky A.S., Chaos Solit. Fract. 5, 1893 

(1995) . 

[14] Boccaletti S., Bragard J., Arecchi F.T., Phys. Rev. E. 
59, 6574 (1999). 

[15] Bragard J., Arecchi F.T., Boccaletti S., Int. J. Bifurca- 
tion and Chaos 10, 2381 (2000). 

[16] Ticos CM., Rosa E., Pardo W.B., Walkenstein J.A., 
Monti M., Phys.Rev.Lett. 85 (2000). 



[17] Ahlers V., Pikovsky A., Phys.Rev.Lett. 88, 254101 
(2002). 

[18] Valladares D.L., Boccaletti S., Feudel F., Kurths J., 

Phys. Rev. E 65, 055208 (2002). 
[19] Mendoza C, Boccaletti S., Politi A., Phys. Rev. E. 69, 

047202 (2004). 

[20] Pecora L.M., Carroll T.L., Jonson G.A., Mar D.J., Chaos 
7, 520 (1997). 

[21] Boccaletti S., Valladares D.L., Phys. Rev. E 62, 7497 
(2000). 

[22] Anishchenko V.S., Vadivasova T.E., Journal of Commu- 
nications Technology and Electronics 47, 117 (2002). 

[23] Anishchenko V.S., Vadivasova T.E., Journal of Commu- 
nications Technology and Electronics 49 (2004). 

[24] Hramov A.E., Koronovskii A. A., Chaos 14, 603 (2004). 

[25] Osipov G.V., Hu B., Zhou C, Ivanchenko M.V., Kurths 
J., Phys.Rev.Lett. 91, 024101 (2003). 

[26] Romano M.C., Thiel M., Kurths J., von Bloh W., Phys. 
Lett. A 330, 214 (2004). 

[27] Koronovskii A.A., Hramov A.E., JETP Lett. 79, 316 
(2004). 

[28] Daubechies I., Ten lectures on wavelets (SIAM, 1992). 
[29] Kaiser O, A friendly guide to wavelets (Springer Verlag, 
1994). 

[30] Koronovskii A. A., Hramov A.E., Continuous wavelet 
analysis and its applications (In Russian) (Moscow, Fiz- 
matlit, 2003). 

[31] Grossman A. and Morlet J., SIAM J. Math. Anal. 15, 
273 (1984). 

[32] Soohoo R.F., Microwave Electronics (Addison- Wesley 

Longman, 1971). 
[33] Infrared and Millimeter Waves (Academic Press, New 



11 



York, 1983), K. Button ed. 
[34] Granatstein V.L., Alexeff I., High Power Microwave 

Sources (Artech House Microwave Library, 1987). 
[35] Felch K.L., Danly B.G., Jory H.R. et al, Proceedings 

IEEE 87, 752 (1999). 
[36] Bezruchko B.P., Kuznetsov S.P., Trubctskov D.I., JETP 

Lett. 29, 162 (1979). 
[37] Ginzburg N.S., Kuznetsov S.P., Fedoseeva T.N., Radio- 

phys. Quantum Electron 21, 728 (1979). 
[38] Levush B., Antonsen T.M., Bromborsky A., Lou W.R., 

Carmel Y., IEEE Trans. Plasma Sci. 20, 263 (1992). 
[39] Nusinovich G.S., Vlasov A.N., Antonsen T.M., 

Phys.Rev.Lett. 87, 218301 (2001). 
[40] Dronov V., Hendrey M.R., Antonsen T.M., Ott E., Chaos 

14, 30 (2004). 



[41] Trubetskov D.I., Mchedlova E.S., Anfinogentov V.G., 
Ponomorenko V.I., Ryskin N.M., Chaos 6, 358 (1996). 

[42] Trubetskov A.E., Hramov A.E., Lectures on microwave 
electronics for physicists (In Russian) (Vol. 1. Fizmatlit, 
Moscow, 2003). 

[43] Packard N.H., Crutchfield J. P., Farmer J.D., Shaw R.S., 

Phys. Rev. Lett. 45, 712 (1980). 
[44] Takens F., in Lectures Notes in Mathematics, edited by 

Rand D., Young L.-S (NY.: Springler-Verlag, 1981), p. 

366. 

[45] Schwarz U., Kurths J., Kliem B., Kriiger A., Urpo S., 
Astron. Astrophys., Suppl. Ser. 127 (1998). 

[46] Wavelets in Physics (Cambridge University Press, 2004), 
J.C. Van den Berg ed. 



